Imaging ambipolar diffusion of photocarriers in GaAs thin films 
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Images of the steady-state luminescence of passivated GaAs self-standing films under excitation by 
a tight ly-fo cussed laser are analyzed as a function of light excitation power. While unipolar diffusion 
of photoelectrons is dominant at very low light excitation power, an increased power results in a 
decrease of the diffusion constant near the center of the image due to the onset of ambipolar diffusion. 
The results are in agreement with a numerical solution of the diffusion equations and with a physical 
analysis of the luminescence intensity at the centre of the image, which permits the determination 
of the ambipolar diffusion constant as a function of electron concentration. 

PACS numbers: 



I. INTRODUCTION 

Ambipolar diffusion is the term used to describe the 
diffusion of electrons and holes in semiconductors when 
their respective concentrations are such that the elec- 
trostatic coupling between the two populations can no 
longer be neglected. From a practical viewpoint, this 
phenomenon must be accounted for when designing any 
bipolar device. After the initial work on electrostatic 
coupling between electrons and holes- 1 , significant the- 
oretical and experimental work has been published on 
ambipolar diffusion in bulk materials 2,3 as well as in 
heterostructures^£. The majority of recent studies con- 
sider undoped material so that the ambipolar diffusion 
constant is only related to hole diffusion 6 8 or to exci- 
tonic transport^. The dependence of the ambipolar dif- 
fusion constant, D a = (D n a p + D p a n )/(a p + a n ), on the 
unipolar diffusion constant D n (D p ) of electrons (holes) 
and of their partial conductivities (J n ((J p ) has never been 
detailed experimentally. Furthermore, the effect of the 
electric field induced by spatial separation of electrons 
and holes has never been evaluated precisely. 

Here we present an optical investigation of ambipo- 
lar diffusion of photoexcited carriers in a thin slab of p + 
GaAs (3 fim thickness) passivated on both sides by 50 nm 
thick GalnP layers (see Fig. [1]). The sample is excited at 
its center by a tightly-focused laser along the z direction 
such that steady-state imaging of the luminescence in- 
tensity enables us to monitor the diffusion profile of mi- 
nority carriers^. The resulting profiles are interpreted 
using two distinct and complementary approaches: i) a 
numerical resolution of the coupled diffusion equations 
for electrons and holes, and ii) a simple qualitative es- 
timate of the electron concentration at the center which 
yields the power dependence of the luminescence thereby 
permitting D a to be evaluated as a function of photo- 
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FIG. 1: Principle of the experiment: A thin, self-supported 
(3/im) GaAs sample is excited by tightly- focused above 
bandgap light (red arrows and top, left inset). An image of 
the bandgap emission is monitored (blue arrows and top, right 
inset). Since the surface recombination is quenched by thin 
GalnP films above and below the GaAs, this image reveals 
the diffusion of carriers within the GaAs. 



electron concentration. 



II. AMBIPOLAR DIFFUSION OF CARRIERS 
IN A THIN SEMICONDUCTING SLAB 



A. Coupled diffusion equations 



In photo-excited p + GaAs the drift-diffusion equations 
for electrons and holes are 
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g - K(N A + 5p)n + V ■ fi n nE + D n Vn (1) 
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: g-K(N A + Sp)n+V- -fi p (N A + Sp)E 



D p VSp 

where Sp is the concentration of photogenerated holes 
and N A is the concentration of acceptors which (in the 
following discussion) will be assumed to be fully ionized. 
K is the bimolecular electron-hole recombination coef- 
ficient and fi n and ji p are the electron and hole mobil- 
ities respectively. Non-radiative bulk recombination is 
neglected for the purposes of this discussion. The terms 
involving the electric field (E) in Eqs. [T] and [2] are respon- 
sible for the electrostatic coupling between electrons and 
holes. In this case E is the internal electric field resulting 
from the spatial distribution of electrons and holes. It is 
given by the Poisson equation 
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where e is the permittivity and q the absolute electronic 
charge. By equating Eqs. Q] and [2] in steady-state, an in- 
dependent expression for the electric field in terms of the 
diffusion constants and concentration gradients can be 
obtained. Using this one may re-write the drift-diffusion 
equation for electrons in the form: 

= g - K(N A + 5p)n + V • Lo a Vn - D' a V(n - Sp)} (4) 



where 



D n Vp(N A + Sp) + D P n n n 
Unn + fi P (N A + Sp) 
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is the usual value of the ambipolar diffusion constant and 
D' a = D p jjL n n/(iJL n n + /jl p (N a + Sp)) gives the magnitude 
of the correction due to the local departure (Sp — n) from 
charge neutrality. The spatial distributions of electron 
and hole concentrations are finally calculated using Eq. 
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Since ambipolar diffusion will be evaluated by varying 
the incident light power and hence the photoelectron con- 
centration, the effect of Fermi blockade on the diffusion 
constants should also be accounted for. In this case the 
diffusion constant depends on the electron concentration 
via the position of the quasi- Fermi level, Ep e , when the 
photo-electron concentration becomes comparable with 
the effective density of states of the conduction band (i.e. 
when the electron gas becomes weakly degenerate). The 
electron diffusion constant is then written 
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where F n ((j)) = J °° x n (exp(x — (j)) + l) _1 dx and = 
jiinksT/q is the low concentration (non-degenerate) value 
of the electron diffusion constant. Here k B is the Boltz- 
mann constant and T is the temperature. The Fermi 
energy is related to the electron concentration in the con- 
duction band by n = x n p(<p)(exp((j) — Ep e ) + l) -1 d</> 
where p((j>) is the density of states in the conduction band 
at energy (j). 
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FIG. 2: (a) The normalized luminescence cross section for a 
light excitation power of 0.013 mW, 0.096 mW, 0.24 mW, 0.49 
mW, and 1.4 mW (curves a to e respectively). Curve f is the 
laser profile, (b) Self consistent calculations under the same 
conditions. In both cases the shape of the profile at low power 
and large r yields a diffusion length of 21.5 /mi. The effect 
of an increase of light intensity is to decrease the diffusion 
length near r = due to ambipolar diffusion, while the slope 
of the logarithmic plot remains practically unchanged at large 
r where u — u/Na is small, (c) Calculated spatial dependence 
of the normalized electronic concentration. The top frame 
shows the low power case (curve a of the right panel), while 
the bottom frame shows the high power case (curve e of the 
right panel). 



B. Power dependence of the luminescence at the 
center. 



It is shown here that simple estimates of the electronic 
concentration at the center of the image (r = 0) can be 
used to qualitatively investigate the unipolar and am- 
bipolar diffusion regimes. For relatively low powers it is 
assumed that the effect of degeneracy on the Einstein re- 
lation is weak, so that the low concentration value of 
the diffusion constant can be used. As will be verified a 
posteriori, at the center (r = 0) it is reasonable to assume 
charge neutrality (n = Sp) so that Eq. [H only contains 
the generation, recombination and ambipolar diffusion 
terms. Secondly, since the diffusion length is larger than 
the Gaussian width a of the laser spot, diffusion domi- 
nates bulk and surface recombination and is thus the de- 
termining factor for the steady-state photoelectron con- 
centration at the center. Assuming that diffusion parallel 
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and power are u = u/Na and 

r)P(l - R) 
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Here P* = irhudD°N A /[r](l - R)} is a power. The solu- 
tion of this equation is 



Au = p(l + 0) - 1 + \/8p+[l-p(l + /3)] 2 



(11) 



As seen from the above approximations, and using Eq. [TJ 
the luminescence at the center is proportional to u(l+u). 




Incident Power ( mW) 

FIG. 3: Luminescence at r = 0, normalized with respect to 
the light excitation power as a function of incident power (red 
squares). Curve a is the adjustment calculated using Eq. [TT] 
and curve b is the result of the ab-initio calculation. 



to the surface can be characterized by a rate r^, one has 
a 2 = DaTdC 1 = D aT * d (7) 

where £ is a numerical factor close to unity. At low power, 
using the value of the unipolar diffusion constant, one 
finds « 6 x 10 -11 s, i.e. about three orders of magni- 
tude smaller than the typical photoelectron lifetime r of 
p+ GaAsii. 

After generation over a characteristic depth 1/a, where 
a is the absorption coefficient, the photoelectrons and 
holes undergo lateral diffusion as well as diffusion along 
z. Since the latter does not change the luminescence 
intensity, the one dimensional diffusion along z can be 
treated independently of the lateral, two dimensional dif- 
fusion. This permits the computation of quantities av- 
eraged along z over the thickness d of the sample. The 
average rate of creation of the photoelectron concentra- 
tion is given by 



g = P(l- R)(/(his7ra 2 d) = Cg* 



(8) 



where R is the reflectivity of the sample surface and ( 
is a numerical factor close to unity. Considering n to 
be homogeneous as a function z, its value is given by 
n = T]g*T% where r] = summarizes the above approxi- 
mations. Using Eqs. [5j[7]and[8j one finds that n does not 
depend on the size of the laser spot a. It is the solution 
of the second degree equation 



2^ 2 + [l-p(l + /3)]ii-p = 
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where f3 = ji n / fi p and the reduced values of concentration 



III. EXPERIMENTAL 

The samples are p + beryllium-doped GaAs thin films 
(Na ~ 10 1T cm -3 ) of thickness 3 jam assembled onto SiC 
substrates^. Reduction of the surface recombination ve- 
locity is provided by 50 nm thick layers of Gao.51Ino.49P 
deposited on each face of the film. The samples are ex- 
cited by a tightly- focussed laser beam of energy 1.59 eV 
in a modified Nikon Optiphot 70 microscope. The laser 
profile, shown in curve f of Fig. [TJ is close to a Gaussian 
profile exp[— r 2 /<j 2 ] with a w 0.93/im. The luminescence 
cross sections were recorded using an appropriate filter 
in order to filter out the excitation wavelength. Curve 
a was taken using a very low excitation power (13 fiW). 
The spatial extent of these profiles is much larger than 
that of the laser, thus revealing electron diffusion in the 
film. At this very low power, diffusion is assumed to 
be unipolar, and the whole cross section is interpreted 
using a diffusion length L = 21.3/im 13 , and is nearly ex- 
ponential for r > 12/im. Shown in curves b, c, d and e 
of Fig. [2] are the spatial dependences of the normalized 
cross sections for increasingly high powers up to 1.4 mW, 
above which the luminescence spectrum reveals a heat- 
ing of both the electron gas and the lattice. Curves b, c, 
d, and to some extent e show little change of the profile 
slope far from the center where the photoelectron con- 
centration is small, thus revealing that the electron dif- 
fusion constant at large radii is close to its unipolar value. 
On the other hand at small radii, the slope strongly in- 
creases indicating a reduction of diffusion constant to its 
ambipolar value given by Eq. [5] for large photoelectron 
concentrations. 

Shown in Fig. [3] is the luminescence magnitude at the 
center, /pl(0) normalized to the incident power and to 
a value of 1 at the lowest power value. This signal is 
close to unity up to about 0.1 mW and reaches values 
larger than 20 for the maximum excitation power. The 
relative excess of carriers at the center is consistent with 
the decrease of diffusion constant due to the progressive 
onset of ambipolar diffusion. 



IV. INTERPRETATION 

In order to interpret the above results, the photoelec- 
tron lifetime r was first measured using time-resolved mi- 



4 



crowave conductivity 14 and found equal to 30.7 ns. The 
good correspondence with the radiative recombination 
time for the nominal doping leveU 3 - is further proof that 
nonradiative surface and bulk recombination processes 
are negligible. Since L is known at very low excitation 
power (curve a of Fig. 2), the value of the unipolar dif- 
fusion constant can be estimated to be « 150 cm 2 /s. 
Finally, as found from the literature 11 one has /3 « 10 so 
that the relevant quantities describing charge diffusion 
are known. 



A. Electronic concentration and luminescence 
intensity at r — 0. 

The normalized luminescence intensity shown in Fig. 
[3] is given by 



WO) 
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where u is the normalized electronic concentration de- 
fined in Sec. II, and Po is the smallest experimental 
power value, corresponding to u = uq. Shown in Fig. [3] is 
the calculated power dependence of /pl(0), using P* =2 
mW. Very good agreement is then obtained using Eq. 
[TTJ Na = 1 x 10 17 cm -3 and a power-independent value, 
close to unity, of rj ~ 1.25. This justifies the main physi- 
cal, but not completely trivial, approximations made for 
obtaining the expression for n. 

The calculated power dependences of the reduced val- 
ues of the ambipolar diffusion constant D a /D^ and of 
the luminescence intensity u(l + u) are shown versus u in 
curves a, c, and d of Fig. HJ Switching from the unipolar 
to the ambipolar regime is revealed by the decrease in 
the diffusion constant. For the maximum value of u one 
finds D a /D^ ^0.2. This result is in agreement with Eq. 
[5j which gives D a « 2D p « 2f3~ 1 D^ L in the limit where 
n ^> Na- It is also seen that u increases faster than the 
light power and that its value at maximum power is of 
the order of 3Na- The power dependence of the lumines- 
cence intensity starts to differ from that of the electron 
concentration for P « 0.1 mW. 

Since the electron concentration at high power is com- 
parable with the intrinsic density of states in the con- 
duction band, the effect of the concentration dependence 
of the electron diffusion constant, described by Eq. [6j 
needs to be evaluated. To first order, taking D n of the 
form D n = D^(l + n/n ), one finds n ~ 1.2 x 10 18 
cm -3 . Eq. [9] becomes a third degree equation includ- 
ing the parameter Na/tlq. Shown in curve b of Fig. [4] 
is the resulting power dependence of the ambipolar dif- 
fusion constant. D n only differs from (curve a) for 
powers larger than about 0.3 mW. For the maximum 
power (P = 1.4 mW) the increase in D n gives a value 
of u slightly smaller than that shown in curve c and cor- 
responds to (D n - P>°)/P>° w 12 %. Given that this 
marginal increase is not unambiguously evident from the 



data, it is reasonable to take a concentration-independent 
electronic diffusion constant. 



B. Luminescence profiles 

In a separate, complementary approach, the coupled 
equations, Eq. [2j Eq. [3] and Eq. 2] were solved self- 
consistent ly using a commercial finite element package. 
This yields the electronic concentration and the photolu- 
minescence intensity at all positions within the sample. 
The bottom panel of Fig. [2]shows the normalized maps of 
electronic concentrations near the center for the smallest 
and for the largest power. It is first verified that at r = 
the relative variation of n as a function z is of the order 
of 40% at small power and of 50% at large power. This a 
posteriori justifies the assumption of homogeneous con- 
centration as a function of depth taken Sec. IIB. Further- 
more since L ^> the variation of n(r, z) as a function of 
z is quite weak in both cases as soon as r is comparable 
with the thickness (note that the horizontal scale in Fig. 
[2^ is much smaller than in Fig. [2ji). For numerical calcu- 
lation of the photoluminescence profiles as a function of 
r it is therefore not a bad approximation to take z = d/2. 

The calculated luminescence profiles are shown in the 
right panel of Fig. [2] after normalization at r = for 
Na — 2.5 x 10 1T cm -3 . The overall behavior of the ex- 
perimental profiles is correctly interpreted by the model 
described above, although slight differences between the 
ab-initio calculations and the experimental results are 
apparent. This is most evident at high power where the 
shape of the profile depends very sensitively on the re- 
duced concentration, u (i.e. on the exact doping density 
and on the incident power). Any small variation in Na 
(whose exact value is not known) or in the incident power 
results in a large relative variation of the luminescence in- 
tensity for large r. For example, the use of Na = 1 x 10 17 
yields an r dependence of the normalized luminescence 
profile that is far too strong. Curve b of Fig. [3] shows 
^pl(O). In the case of the ab-initio calculation, the ratio 
is calculated after integration over the whole thickness of 
the sample and over a lateral radius of the order of that 
of the excitation spot. Once again, although the quali- 
tative shapes of the calculated and experimental curves 
are in reasonable agreement, there are quantitative differ- 
ences between the curves. As above, this is particularly 
so at high power where the luminescence intensity de- 
pends sensitively on u. Undoubtedly better agreement 
could be obtained by varying several parameter values 
(Na, ft, etc. . . ) but doing so is tedious and not partic- 
ularly revealing from a physical point of view. It is also 
possible that the slight difference is due to photon recy- 
cling which could yield a luminescence profile somewhat 
larger than that due to carrier diffusion alon o 15 i 16 . 

One advantage of the ab-initio calculation is that it 
can be used to evaluate the assumption of local charge 
neutrality (i.e. n = Sp) that is made in all discussions of 
ambipolar transport 1 . Fig. [5] shows the spatial distribu- 
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FIG. 4: The calculated dependence of the diffusion coefficient 
normalized to its unipolar value as a function of light excita- 
tion power. Curve a(b) corresponds to the result calculated 
without (with) Fermi blockade taken into account. Also shown 
are to dependences of the electronic concentration u normal- 
ized to the acceptor concentration (curve c) and of u(l + u) 
(curve d) which is proportional to the luminescence intensity. 
The results of the full numerical calculation (green squares) 
are in excellent agreement with the simplified analytical ap- 
proach. Shown for comparison is a straight, dotted line of 
slope unity. 



the excellent agreement obtained between the exact nu- 
merical and approximate analytic calculations of D a /D^ 
at r — shown in Fig. |4j 
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FIG. 5: Estimation of the internal electrical field and the va- 
lidity of the local neutrality approximation for an incident 
powers of 0.013 mW (curves a' and b') and 1.4 mW (curves a 
and b). Curves a and a' show the spatial distribution of the 
relative difference between electron and hole concentration. 
Curves b and b' show the relative difference between the spa- 
tial distribution of electrons obtained when n — n — Sp in 
Eq. [Hand when the electric field is accounted for. 



tion of the relative difference (n — Sp)/n at z = d/2 for the 
lowest (curve a') and highest (curve a) incident powers. 
In both cases there is an excess of holes near r = and 
a compensating excess of electrons at a distance larger 
than 3-4 /am. As expected, the relative excess of holes at 
the center is larger at low powers where ambipolar dif- 
fusion is absent. In the presence of ambipolar diffusion, 
electrons and holes have a tendency to diffuse together 
and the relative difference drops by a factor of 10. Since 
the permittivity e in Eq. [3] is small, these observations 
do not necessarily imply that the term proportional to 
E is negligible. In order to validate the assumption of 
local charge neutrality, the electronic concentration n' 
obtained when neglecting the last term of Eq. [4] is calcu- 
lated. Shown in curves b and b' of Fig. [5] is the relative 
value (n' — n)/n for the highest and lowest powers respec- 
tively. Unsurprisingly, the term proportional to E is more 
important at higher power , although at worst, assuming 
n — Sp introduces an error of the order of 10% into the 
resulting concentration profiles. More importantly, the 
error is smallest at r = 0, indicating that the simplifying 
assumptions used above to analyze the luminescence in- 
tensity at the center are reasonable. This is confirmed by 



V. CONCLUSION 

Imaging of the luminescence profile created by a highly 
focused excitation and emitted by a 3/im thick p + 
GaAs clearly reveals ambipolar diffusion as the excita- 
tion power is increased. The switching from unipolar 
to ambipolar diffusion of photocarriers is investigated as 
a function of electron concentration and the results are 
analyzed using a numerical resolution of the coupled elec- 
tron and hole diffusion equations, as well as the Poisson 
equation. It is found that the effect of the electric field 
induced by ambipolar diffusion can be significant away 
from the center. In contrast, this effect is reduced near 
the center so that a simple calculation of the power de- 
pendence of the luminescence intensity can be performed. 
The results are interpreted using a single parameter, de- 
fined in Eq. [10] as a power P*, which depends on acceptor 
concentration, slab thickness and unipolar electron diffu- 
sion constant. The experimental results at the center 
are in very good agreement with the predictions of this 
model, using a reasonable value of P*. 
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